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Abstract 

We study the contribution of gluon induced partonic subprocesses to Z'j pair pro¬ 
duction at hadron colliders. These processes contribute only at next-to-next-to-leading 
order but are potentially enhanced by two factors of the gluon parton densities. How¬ 
ever, we hnd that their contribution is modest and that next-to-leading order calcula¬ 
tions give reliable predictions. 


1 Introduction 


The substantially increased number of Z 7 pairs that are expected to be produced 
in Run II at the Tevatron will allow an improvement to previous studies i i i of 
Zy-pair production at hadron colliders. Accordingly, the ZZ^ and Z 77 couplings, 
which are absent in the Standard Model, can be investigated in more detail. Assuming 
no deviations from the Standard Model are found, tighter limits on such anomalous 
couplings will be possible. 

In order to fully exploit the experimental data it is important to have theoretical 
predictions of sufficient accuracy. Next-to-leading order cross sections for this process 
have been calculated quite some time ago [^. Subsequently, the theoretical predic¬ 
tions have been improved by taking into account the leptonic decay of the Z-boson, 
anomalous couplings and full leptonic correlations iiBi- 

For some distributions, the next-to-leading order corrections have been found to be 
very large. In particular in the large transverse momentum region, which is interest¬ 
ing for investigating anomalous couplings, the corrections can increase the theoretical 
prediction of the cross section by several 100%. The reason for this large corrections 
is that at NLO new partonic processes contribute. Whereas the tree-level calculation 
only incorporates the partonic process qq Zj —> iij, at next-to-leading order also 
processes with a gluon in the initial state such as qg —>• Zyg —> Hq^j contribute. The 
large gluon density at low x overcomes the suppression by an additional factor of 
Thus, even though this contribution is formally next-to-leading order, numerically it 
can be as important as the leading-order term. 

This immediately leads to the question whether at NNLO a similar effect occurs. 
Higher-order corrections to partonic processes that are already present at NLO are 
not expected to lead to large corrections. However, at NNLO a new class of partonic 
processes have to be taken into account, namely processes with two gluons in the initial 
state. Such processes are suppressed by two factors of a* but are potentially enhanced 
by two factors of the gluon parton distribution function. If there is a similar effect as 
at NLO, this could lead to corrections as important as the Born term. 

This question has been studied for WZ and W 7 production 


10 . For these 


processes it has been found that the gluon induced NNLO corrections are negligible. 
However, there is an important difference in Z 7 production with respect to WZ, W 7 
production. In fact, there are generally two classes of contributions for vector-boson 
pair production with two gluons in the initial state. There are tree-level processes 
gg VVqq and loop processes gg —> VV, where V denotes some vector boson. In 
the case of hFZ and hFy the loop processes vanish due to charge conservation but for 
Z 7 this is not the case. The study in |P, showed that the tree-level processes result 
in a very small contribution for WZ and W'y. While it is natural to expect that this 
is also the case for Z 7 production, for this process there remains the question of how 
important the gluon induced loop corrections are. 

The loop induced process gg —> Z 7 has beeen studied previously [0, |T^ |T^, where 
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it has been found that these contributions can be quite large but are not dominant. 

This work extends the previous analyses in that both classes of gluon induced 
corrections to Z 7 production are included and anomalous ZZ^ and Z 77 couplings as 
well as full leptonic correlations are taken into account. Even though the gluon induced 
correction is substantially more important for Z 7 production than for hhy, WZ \t can 
still be safely neglected. In fact, not only are the contributions from tree-level processes 
tiny (as in the IE 7 , WZ case) but also the loop processes are smaller than anticipated. 
The latter effect is entirely due to a change in the parton distribution function used in 
our analysis compared to 
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2 Calculation 


Our aim is to compute the part of the NNLO corrections to Z 7 pair production that 
is enhanced by two factors of the gluon parton distribution function. This consists of 
two parts. Firstly, we have to compute the tree-level amplitudes gg —> qqZ^ —> qqli'y 
and integrate over the quark-antiquark hnal state phase space. Secondly, we have to 
evaluate the loop diagrams gg ^ Z'j —>■ ii'j. 

The contributions of the tree-level diagrams gg —>■ qqZj —>■ ggHj can be com¬ 
puted in a very similar manner as for gg —> gglV-y, presented in |^. This is a double 
bremsstrahlung NNLO contribution and in general is difficult to integrate over the 
phase space. However, due to the simple structure of the soft and collinear limits of 
the corresponding matrix elements, the integration over the singular regions of the 
phase space is much simpler than for a general NNLO calculation. In fact, as dis¬ 
cussed in [ 0 , the only singularities we are concerned with are single and double initial 
state collinear singularities. They come from the region of phase space where the two 
incoming gluons independently split into a quark-antiquark pair. 

In order to perform the integration over phase space we use a generalization of 
the subtraction method presented in The singularities that arise upon integration 


over the region of phase space with the hnal state (anti)quark collinear to the incoming 
gluons are absorbed into the parton distribution functions. It should be noted that this 
contribution is factorization scheme dependent. Since we use parton densities obtained 
in the MS'-scheme we evaluate the partonic cross section in this scheme. We should 
also mention that the factorization scheme dependence is beyond the approximation 
we make since it is not enhanced by two factors of the gluon density. 

The tree-level diagrams for gg qqVV have been computed previously [|T^. We 
included anomalous couplings ZZ'j and 2 ^ 77 . They are incorporated by using the 
Za{qihg{q 2 )Zi,{p) vertex [T|, 0 


^%z{Qi,Q2,p) = 


z{p^ - qj) 
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hf{q^g^^-q^g>^^) + 


ho 
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where it should be noted that the Z-boson with momentum p is off-shell and the other 
two gauge bosons are on-shell. The vertex for Za{qi)'Ji3{q2)'Jn{p) is as in eq. (|I[), but 
with —>■ 0 and hf ^ h]. 

The inclusion of anomalous couplings also necessitates the recalculation of the gluon 
induced loop diagrams. First of all, there are box diagrams. They have been calculated 
some time ago BBig but these results were obtained summing over helicities and 
did not include the decay of the Z-boson. However, since the axial coupling of the Z 
does not contribute, the amplitude can also be extracted from (the Np-part of) the 
QCD amplitudes gg ggy |^, [T^ . In fact we simply have to change the gluon that 
decays into the qq pair into a Z. This can be done by picking the relevant subamplitudes 
of the gg ggy process and changing the couplings accordingly. 

In addition to the box diagrams there are triangle diagrams with a triple ZZ'y or 
Z'j'j vertex. These diagrams are only present for anomalous couplings. Furthermore, 
due to Yang’s theorem this contribution is only non-vanishing because of off-shell 
effects. It should also be stressed, that the vertex of eq. (0) does not include all 
possible anomalous couplings. The vertex in eq. m does not include terms ~ p^ 
since it is assumed that the gauge bosons couple to a conserved current. However, in 
our situation this is not the case and aditional anomalous couplings are possible . 
There is one additional CP conserving {p^qiq^^avap) and two additional CP violating 
couplings {Pfiga/s and —g2)a(p —<?i)/3)- In our calculation we have not implemented 
these additional anomalous couplings, since — given the small size of the gluon induced 
corrections — it is very unlikely that useful information about these couplings can be 
extracted. 

The effects of the top quark are taken into account by performing an expansion in 
l/rrit and neglecting terms of order This approach is motivated by assuming 

that the partonic center of mass energy y/§ < rrit. Even though this is not generally 
true at the LHC, the fact that the gluon distribution is strongly peaked towards small 
values of x, makes this a reasonably good approximation. 

To deal with the problem of photon isolation, we will use the procedure introduced 
by Frixione |^, i.e. we reject all events unless the transverse hadronic momentum 


deposited in a cone of size Rq around the momentum of the photon fulhlls the following 
condition 

Y.PT.e{R- < PT ^), ( 2 ) 

i vl - cosR qJ 

for all R < Rq, where the ‘distance’ in pseudorapidity and azimuthal angle is dehned 
by Ri^ = -F ((/>* - 0^)2. 


3 Numerical Results 


For the numerical results presented in this section we follow closely the analysis pre¬ 
sented in p]. We use the MRST 2001 parton distribution functions with the 
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one-loop expression for the coupling constant {as{Mz) = 0.119). Strictly speaking, 
for the gg induced processes we should use a NNLO parton distribution. However, 
it is not expected that a change to the not yet fully available NNLO distributions 
would alter o ur conclusio ns. The factorization and renormalization scales are hxed to 
hF = hi? = ^M| + (p^)2. 

The mass of the Z has been set to Mz = 91.187 GeV. For the couplings of the 
vector bosons with the quarks we use a = a{Mz) = 1/128 whereas for the photon 
coupling we use a = 1/137. Note that we do not include the branching ratios for the 
decay of the Z into leptons. 

For our numerical results we use a set of standard cuts. Charged leptons are required 
to have pt > 20 GeV and p < 2.5. The photon transverse momentum cut we use is 
pj > 20 GeV, while for the isolation prescription in eq. (@) we set i?o = 1- 



pj |GeV| 

Figure 1: The contribution to the distribution in Zy production at the LHC separated 
into the various partonic initial states: qq (solid line), gq (dashed) and gg. The gluon 
induced contribution is separated into the loop and bremsstrahlung part. 

Since the potential enhancement due to the gluon densities is larger for increasing 
center of mass energy, we restrict ourselves to studying the situation at the LHC. Thus, 
we consider Zy production in proton-proton collision at an energy y/s = 14 TeV. 

The canonical quantity we use to investigate the importance of the gluon induced 
corrections is the transverse momentum of the photon. In Figure we show the con¬ 
tribution of the various initial partonic states to dcx/dpj. The solid line shows the 
contribution at NLO of the initial state qq. This is the dominant part and differs only 
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slightly from the tree-level result. The gq (and gq) initial states enter at NLO and can 
result in a correction as big as 70% which is shown as the dashed line. Finally, there are 
the gg induced processes. We split the result into the loop part (short dashed) and the 
tree-level bremsstrahlung part (dotted). The latter is negative and very small. This is 
very similar to what has been found for WZ and fFy production 0. The loop part 
is substantially bigger but still does not exceed 5%. Note that differs quite a bit from 
what has been found in previous analyses |^, |T^. We checked that this is entirely due 
to an update in the gluon distribution function. In fact, our results agree with those 
in |T^, [T^ if we use the same parton distribution functions. 

The situation does not change if anomalous couplings are added. The relative 
importance of the gluon induced processes is still small. Thus, there is no hope of 
gaining any information about non-standard anomalous couplings that in principle 
contribute to this process. 

We also investigated other quantities and the picture remains the same. The loop 
corrections have more or less the expected size but are never very important. The 
reason for the smallness of the bremsstrahlung corrections is the following: in the 
important region of small partonic center of mass energy a/^, corresponding to small 
X, the partonic cross section dd is small and turns negative. Only for increasing \/5 we 
hnd that da is of the expected size, i.e. suppressed by two orders of However, in 
this region x is not small anymore and, thus, there is no enhancement due to the gluon 
distribution. This might well be related to the very simple structure of the singularities 
of the partonic cross section. For a process with a more complicated structure, e.g. 
involving t-channel gluon exchange, the bremsstrahlung corrections could easily be 
signihcantly higher. 


4 Conclusions 


We studied the contribution of the partonic processes gg qqZ'j and gg Z'j to 
the production of Zj pairs at hadron colliders. These contributions are enhanced by 
two factors of the large gluon density. This could potentially overcome the 0{ag) 
suppression. However, we found that under no circumstances are these contributions 
particularly important. In fact, they are even substantially smaller than anticipated 
from previous analyses |]T^ O, a change that is due to using updated parton distribu¬ 


tion functions. 

In summary, the next-to-leading order cross section for Z'j production provide us 
with a reasonably precise theoretical prediction. The large NLO corrections due to the 
opening of a new partonic channel with qg initial states are taken into account. Higher 
order corrections to these partonic processes are expected to be well under control and 
the new partonic channels that open at NNLO, namely gg initial state processes, do 
not result in large corrections. 
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